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ABSTRACT 

Using the Sloan Digital Sky Survey (SDSS) galaxy group catalogue of Yang et al. 
(2007), we study the average colour (representing star formation history) and average 
concentration (representing mass assembly history) of satellite galaxies as function of 
(i) their stellar mass, (ii) their group mass, and (iii) their group-centric radius. We 
find that the colours and concentrations of satellite galaxies are (almost) completely 
determined by their stellar mass. In particular, at fixed stellar mass, the average 
colours and concentrations of satellite galaxies are independent of either halo mass or 
halo-centric radius. We find clear evidence for mass segregation of satellite galaxies 
in haloes of all masses, and argue that this explains why satellites at smaller halo- 
centric radii are somewhat redder and somewhat more concentrated. In addition, the 
weak colour and concentration dependence of satellite galaxies on halo mass is simply a 
reflection of the fact that more massive haloes host, on average, more massive satellites. 
Combining these results with the fact that satellite galaxies are, on average, redder 
and somewhat more concentrated than central galaxies of the same stellar mass, the 
following picture emerges: galaxies become redder and somewhat more concentrated 
once they fall into a bigger halo (i.e., once they become a satellite galaxy). This is a 
clear manifestation of environment dependence. However, there is no indication that 
the magnitude of the transformation (or its timescale) depends on environment; a 
galaxy undergoes a transition when it becomes a satellite, but it does not matter 
whether it becomes a satellite of a small (Milky Way sized) halo, or of a massive 
cluster. We discuss the implication of this 'dearth' of environment dependence for the 
physical processes responsible for transforming satellite galaxies. 

Key v^rords: galaxies: clusters: general - galaxies: haloes - galaxies: evolution - 
galaxies: general - galaxies: statistics - methods: statistical 



1 INTRODUCTION 

In the current paradigm of galaxy formation, it is believed 
that virtually all galaxies initially form as disks due to the 
cooling of gas with non-zero angular momentum in viri- 
alized dark matter haloes. During their subsequent evolu- 
tion, star-forming disk galaxies may be transformed into 
non star-forming early-types via a variety of transformation 
mechanisms. Two disk galaxies of roughly equal mass may 
merge to produce a spheroidal galaxy (Toomre & Toomre 
1972; Negroponte & White 1983), or a disk galaxy may be 
transformed into a spheroid due to the cumulative effect of 
many high speed (impulsive) encounters, called harassment 
(Farouki & Shapiro 1981; Moore et al. 1996). When a small 
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halo is accreted by a larger halo its hot, diffuse gas may 
be stripped thus removing its fuel for future star-formation 
(Larson, Tinsley & Caldwell 1980; Kauffmann, White & 
Guiderdoni 1993; Balogh, Navarro & Morris 2000). Follow- 
ing Balogh & Morris (2000) we refer to this process, that 
results in a fairly gradual decline of the satellite's star for- 
mation rate, as strangulation. When the external pressure 
is sufficiently high, ram-pressure stripping may also remove 
the entire cold gas reservoir of the satellite galaxy (Gunn & 
Gott 1972; Quilis, Moore & Bower 2000), resulting in an ex- 
tremely fast quenching of its star formation. Finally, along 
its orbit a satellite galaxy is subject to tidal forces which 
may cause tidal stripping and heating. 

The efficiencies of all these various transformation 
mechanisms are expected to be strong functions of environ- 
ment. Ram-pressure stripping requires a dense inter-galactic 
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medium and is therefore believed to occur predominantly 
in massive clusters. The same holds for galaxy harassment, 
which requires the presence of many, relatively massive satel- 
lite galaxies in order to be efBcient. Galaxy merging, on the 
other hand, is thought to be suppressed in massive clusters, 
and to preferentially occur in group-sized haloes, while tidal 
stripping and heating should play a role in aJl dark matter 
haloes with little dependence on their viriaJ mass. Finally, 
strangulation can result from both ram-pressure and tides, 
and it is difficult to predict how its efficiency would scale 
with halo mass. The various transformation mechanisms also 
have very different outcomes: while mergers and harassment 
arc thought to transform disk galaxies into spheroids, ram- 
pressure stripping and strangulation only affect the star for- 
mation activity of the galaxy without changing its global 
morphology. Therefore, an assessment of galaxy properties 
as function of their environment sheds light on which of 
these transformation mechanisms operates in what environ- 
ment and with what efficiency. 

Numerous authors have investigated the relation be- 
tween environment and morphology (e.g., Hubble 1936; 
Oemler 1974; Dressier 1980; Postman & Geller 1984; Whit- 
more, Gilmorc & Jones 1993; Hashimoto & Oemler 1999; 
Goto et al. 2003; Mcintosh, Rix & Caldwell 2004; Kuehn 
& Ryden 2005; Blanton et al. 2005b), between environment 
and star formation rate (e.g., Hashimoto et al. 1998; Lewis et 
al. 2002; Dommguez et al. 2002; Gomez et al. 2003; Balogh et 
al. 2004a; Tanaka et al. 2004; Kauffmann ct al. 2004; Kclm, 
Focardi & Sorrentino 2005), and between environment and 
colour (e.g., Tanaka et al. 2004; Balogh et al. 2004b; Hogg et 
al. 2004; Wcinmann ct al. 2006a; Blanton & Bcrlind 2007). 
These and other studies have clearly established that galax- 
ies in denser environments are more massive, redder, more 
concentrated, less gas-rich and have an older stellar pop- 
ulation. Within the paradigm outlined above, the trans- 
formation mechanisms thus seem to preferentially occur in 
denser environments. However, a more detailed interpreta- 
tion of these findings is complicated by the fact that various 
galaxy properties are strongly correlated even at a fixed en- 
vironment. An important outstanding question, therefore, is 
which relationships with environment are truly causal, and 
which are just reflections of other more fundamental corre- 
lations that are not environment dependent. 

With the advent of large redshift surveys, such as the 
Two-Degree Field Galaxy Redshift Survey (2dFGRS; Col- 
less et al. 2001) and the Sloan Digital Sky Survey (SDSS; 
York et al. 2000; Stoughton et al. 2002), several authors 
have started to address these questions using multi-variate 
statistics. Most of these studies agree that morphological 
properties are less strongly related to environment than star- 
formation related properties such as colour and emission-line 
fiux (Kaufi^mann et al. 2004; Blanton et al. 2005b; Quintero 
et al. 2006). This is also in line with recent results indicat- 
ing that the relation between star formation indicators and 
environment is at least partially independent of galaxy mor- 
phology (Balogh et al. 1998; Koopmann & Kenney 1998; 
Poggianti et al. 1999; Christlein & Zabludoff 2005) . 

A physical interpretation of all these results in terms 
of the transformation mechanisms discussed above is ham- 
pered by two shortcomings. First of all, in most studies of 
environment dependence, especially those based on large 
redshift surveys such as the 2dFGRS or the SDSS, one 



does not split the galaxy population into central galaxies 
(those at rest at the centers of their dark matter haloes) 
and satellite galaxies (those orbiting a central galaxy) . This 
is important since most of the transformation mechanisms 
discussed above (tidal stripping and heating, strangulation, 
ram-pressure stripping, and harassment) mainly operate on 
satellite galaxies. However, satellite galaxies only maJte up 
roughly 20 to 40 percent of the entire galaxy population 
(Mandelbaum et al. 2006; van den Bosch et al. 2007a; Tin- 
ker ct al. 2007; Cacciato ct al. 2008). Therefore, unless one 
specifically focuses on satellite galaxies, the environmental 
impact of the aforementioned processes may not be appar- 
ent. 

A second problem that hampers a physical interpre- 
tation of the results is the use of non-intuitive environ- 
ment indicators. Most studies parameterize 'environment' 
through the projected number density of galaxies, either 
within some fixed metric aperture, or out to the distance 
of the n**^ nearest neighbor (with n typically in the range 
5-10). Although these environment indicators are relatively 
easy to measure from large galaxy redshift surveys, they axe 
difficult to interpret in terms of more physical quantities 
(see discussions in Kauffmann et al. 2004 and Weinmann et 
al. 2006a). Arguably the most natural environment indica- 
tors are the halo mass and the halo-centric distance, with the 
virial radius flvir being the natural scale of environment de- 
pendence. This is supported by various observational stud- 
ies which have shown that the environment dependence of 
galaxies seems to be restricted to scales R < i?vir (e.g.. Mo et 
al. 2004; Kauffmann et al. 2004; Quintero et al. 2005; Blan- 
ton ct al. 2006; Blanton & Bcrlind 2007). Galaxy properties 
as function of halo mass can be studied using either clus- 
tering data or using galaxy group catalogues. The former 
is less direct as it requires modelling in order to translate 
an observed correlation function into a description of halo 
occupation statistics. Nevertheless, a great deal of progress 
in this area has been made in recent years (e.g., Yang ct 
al. 2003; van den Bosch et al. 2003, 2007a; Scranton 2003; 
Magliochetti & Porciani 2003; Zehavi et al. 2005; Tinker et 
al. 2005; Colhster & Lahav 2005; Li et al. 2006; Wang et 
al. 2007a; Cacciato ct al. 2008). Galaxy group catalogues 
allow a more direct study of the relation between the prop- 
erties of galaxies and their dark matter haloes. In particular, 
it assigns a halo to each individual galaxy, whereas the clus- 
tering method only assigns haloes to galaxies in a statistical 
sense. Another advantage of using galaxy group catalogues is 
that it is straightforward to separate the galgixy population 
into centrals and satellites and to compute their halo-centric 
distances (in projection). 

In van den Bosch ct al. (2007b, hereafter Paper I), we 
used a large galaxy group catalogue constructed from the 
SDSS to study the impact of satellite specific transforma- 
tion processes. We compared the colours and concentrations 
of satellites galaxies to those of central galaxies of the same 
stellar mass, adopting the hypothesis that the latter are 
the progenitors of the former. On average, satellite galax- 
ies were found to be redder and more concentrated than 
central galaxies of the same stellar mass, indicating that 
satellite specific transformation processes do indeed oper- 
ate. Central-satellite pairs that are matched in both stellar 
mass and colour, however, show no average concentration 
difference, indicating that the transformation mechanisms 
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operating on satellites affect colour more than morphology. 
We also showed that the colour and concentration differ- 
ences of matched central-satellite pairs are completely inde- 
pendent of the mass of the host halo (not to be confused with 
the subhalo) of the satellite galaxy, indicating that satellite 
specific transformation mechanisms are equally efficient in 
host haloes of all masses. In this paper we extend this study 
by investigating the ecology of satellite galaxies. In par- 
ticular, we investigate how the colours and concentrations 
of satellite galaxies arc correlated with their stellar mass, 
their halo mass, and their (projected) halo-centric distance. 
We believe that this is the optimal strategy to investigate 
the effects of satellite-specific transformation mechanisms 
such as strangulation, ram-pressure stripping and harass- 
ment. Throughout we adopt a fiat ACDM cosmology with 
flm = 0.238 and = 0.762 (Spergel et al. 2007) and we 
express units that depend on the Hubble constant in terms 
oih = Ho/wo kms"^ Mpc'^ 



2 DATA 

The analysis presented in this paper is based on the SDSS 
DR4 gala^xy group catalogue of Yang et al. (2007; here- 
after Y07). This group catalogue is constructed applying 
the halo-based group finder of Yang et al. (2005a) to the 
New York University Value-Added Galaxy Catalogue (NYU- 
VAGC; see Blanton ct al. 2005), which is based on SDSS 
DR4 (Adelman-McCarthy et al. 2006). From this catalogue 
we select all galaxies in the Main Galaxy Sample with an ex- 
tinction corrected apparent magnitude brighter than r = 18, 
with redshifts in the range 0.01 < z < 0.20 and with a red- 
shift completeness C > 0.7. As described in Y07, we use this 
sample of galaxies to construct three group samples: Sample 
I, which only uses the 362356 galaxies with measured red- 
shifts from the SDSS, Sample II which also includes 7091 
gala^xies with SDSS photometry but with redshifts taken 
from alternative surveys, and Sample III which includes an 
additional 38672 galaxies that lack a redshift due to fiber- 
collisions, but which we Eissign the redshift of its nearest 
neighbour (cf. Zehavi et al. 2002). Unless specifically stated 
otherwise the present analysis is based on Sample II, which 
consists of 369447 galaxies distributed over 301237 groups. 

The magnitudes and colours of all galaxies are based 
on the standard SDSS petrosian technique (Petrosian 1976; 
Strauss et al. 2002), have been corrected for galactic extinc- 
tion (Schlegel, Finkbeiner & Davis 1998), and have been 
/sT-corrected and evolution corrected (hereafter i^-corrected) 
to z = 0.1, using the method described in (Blanton ct 
al. 2003a). We use the notation '^'^Mr to indicate the re- 
sulting absolute magnitude in the r-band. In addition, for 
each galaxy we compute a stellar mass, M* , using the rela- 
tions between stellar mass-to-light ratio and colour of Bell et 
al. (2003; see Y07 for details). In addition to the magnitude 
and stellar mass, we also use the '^'^{g — r) colour of each 
galaxy as well as its concentration C = rgo / rso . Here rgo and 
rso are the radii that contain 90 and 50 percent of the Pet- 
rosian r-band flux, respectively. In what follows, we interpret 
the '^'^{g — r) colour as an indicator of the galaxy's star for- 
mation history (mean stellar age), though we acknowledge 
that because of dust attenuation this connection may not be 
very tight. In addition, we interpret the concentration pa- 



rameter as reflecting the galaxy's mass-assembly history. As 
shown by Strateva et al. (2001), C is a reasonable proxy for 
Hubble type, with C > 2.6 corresponding to an early-type 
morphology. The link to the mass-assembly history owes to 
the common belief that a galaxy's morphology (in particu- 
lar its disk-to-bulge ratio), is related to the importance of 
major mergers (and other violent disturbances) during its 
mass assembly. 

For each group in our catalogue we have two estimates 
of its dark matter halo mass Mh ; one based on the ranking 
of its total characteristic luminosity, and the other based on 
the ranking of its total characteristic stellar mass. As shown 
in Y07, both halo masses agree very well with each other, 
with an average scatter that decreases from ~ 0.1 dex at the 
low mass end to ~ 0.05 dex at the massive end. In addition, 
detailed tests with mock galaxy redshift catalogues have 
demonstrated that our group masses are more reliable than 
those based on the velocity dispersion of the group members 
(Yang et al. 2005b; Weinmann et al. 2006a; Y07). Further 
support for the accuracy of our group masses comes from an 
analysis of their clustering properties (Wang et al. 2007b). 
In this paper we adopt the group masses based on the stel- 
lar mass ranking. We have verified, though, that none of our 
results change significantly if we adopt the luminosity-rank 
based masses instead (see §4.1) 

Finally, for each group member we determine the pro- 
jected radius -Rproj from the luminosity weighted group cen- 
ter using the angular separation and the angular diameter 
distance at the redshift of the group. We normalize these 
projected radii by the characteristic radius of the group, 
Riso , which is defined as the radius inside of which the dark 
matter halo associated with the group has an average over- 
density of 180. This radius is computed from the mass and 
redshift of the group using equation (5) in Y07. 



2.1 Sample definition 

From Sample II we select all satellite galaxies, defined as 
all group members that are not the most massive group 
member, of groups with an assigned halo mass in the range 
12.0 < log[Mh/lh~^ Mq)] < 15.0. In addition, the satel- 
lites are required to have stellar masses in the range 9.0 < 
\og[M, / {h~^ Mq)] < 11.5 and normalized projected radii 
-Rproj / R 180 < 1.0. Neither of these criteria is particularly 
restrictive, so that the vast majority (~ 88%) of all satel- 
lite galaxies in our group catalogue make it into the sample 
used here. The resulting sample consists of 60245 satellite 
galaxies, which we use to study the interrelations between 
colour, concentration, halo mass, stellar mass, and projected 
halo-centric radius. 

Note that this is not a volume-limited sample. Conse- 
quently, the sample suffers from Malmquist bias, causing an 
artificial increase of the average luminosity (and also stel- 
lar mass) of the satellite galaxies with increasing redshift. 
To correct for this bias, we weight each galaxy by 1/Knax, 
where Vmax is the comoving volume of the Universe out to 
a comoving distance at which the galaxy would still have 
made the selection criteria of our sample. In what follows 
all distributions are weighted by 1/Vmax, unless specifically 
stated otherwise 
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Figure 1. The average satellite mass as a function of the normal- 
ized, projected radius from the group center. Different lines corre- 
spond to different bins in halo mass, as indicated. Two trends are 
evident. First of all, at a given i?proj/Rl80 more massive haloes 
contain, on average, more massive satellites. Secondly, there is 
clear evidence for mass segregation, with more massive satellites 
(in a given halo mass), residing at smaller (projected) radii. 



3 RESULTS 

The main aim of this paper is to investigate the condi- 
tional probability functions P{^'^ {g — r)\Mf , Mu, Rpioj) and 
P(C [M« , Mh, i?proj) of satellite galaxies. In particular, we 
wish to establish which of the conditionals, M«, Mh or 
J?proj/i?i80, is most relevant for setting the colours (related 
to the star formation histories) and concentrations (related 
to the mass-assembly histories) of satellite galaxies. 

In order to facilitate the interpretation of these proba- 
bility functions and their moments we first focus on the in- 
terrelations between the three conditionals. These are nicely 
summarized in Fig. 1, which plots the average stellar mass 
of the satellite galaxies, {log[Af, M©)]), as function 
of 7?proj/i?i80 for different bins in halo mass. Two trends 
are clearly apparent: at fixed /?proj /i^iso , the average satel- 
lite mass increases with increasing halo mass, while at fixed 
Mh, the average satellite mass decreases with increasing 
Rproj/ Rim- The former simply reflects that the charac- 
teristic mass (or luminosity) of satellite galaxies increases 
with increasing halo mass (cf., Yang et al. 2005c; Zheng et 
al. 2005; Skibba, Sheth & Martino 2007; Yang, Mo & van den 
Bosch 2008). The latter indicates that the spatial distribu- 
tion of satellite galaxies in dark matter haloes has undergone 
some mass segregation, and is consistent with the luminos- 
ity segregation observed in galaxy clusters (e.g.. Rood & 
Turnrose 1968; Quintana 1979; den Hartog & Katgert 1996; 
Adami, Biviano & Mazure 1998; Lares, Lambas & Sanchez 
2004; Mcintosh et al. 2005). As we will see below, these two 
trends are essential for understanding the various relations 
between colour and concentration on the one hand and halo 
mass, stellar mass, and halo-centric radius on the other. 



3.1 First Moments 

As a first step in our investigation of Pi^'^{g — 
r)\M*,Mh,Rp,oi) and P{C \M*, Mh, Rproj) of satellite 
galaxies we focus on their first moments. Fig. 2 shows the re- 
lations between the '^'^{g — r) colour of satellite galaxies and 
their halo mass Mh, their stellar mass M*, and their normal- 
ized, projected, halo-centric radius, -Rproj /^?i80 • The upper 
panels show contour plots of the various distributions, while 
the solid lines refiect the average colour. Panel (a) shows that 
the colour distribution of satellite galaxies has a remarkably 
weak dependence on halo mass: over the entire range of halo 
masses probed the colour distribution is clearly skewed, with 
a relatively narrow peak near '^'^{g — r) = 0.85 and an ex- 
tended wing to bluer colours. There is a weak trend, though, 
that the average satellites are somewhat bluer in less massive 
haloes. 

The colour-stellar mass distribution, shown in panel 
(b), reveals a pronounced and narrow red sequence which is 
clearly tilted with more massive red sequence satellites being 
redder. In fact, this colour-stellar mass relation of satellite 
galaxies looks very similar to that for the entire galaxy pop- 
ulation (cf. Fig. 7 in Paper 1). We emphasize that this is 
not a trivial result since satellite galaxies only contribute 
between 20 and 40 percent of the entire galaxy population 
(see §1). As discussed in Paper I, it implies that centrals and 
satellites of the same stellar mass have very similar colour 
distributions, which already puts tight constraints on the 
efficiencies with which the various satellite specific transfor- 
mation processes operate. 

Panel (c) of Fig. 2 shows the distribution of satellite 
galaxies as function of °'^{g — r) and RpToj/Riso- First of all, 
the number of satellites seems to peak around iJproj / Riao = 
0.15; the rapid decline toward smaller radii, however, is an 
artefact of the data sample, and owes to the problem of 
fiber collisions in the SDSS. Since galaxies lost from the sur- 
vey due to fiber collisions do not have any specific colours 
or concentrations, this has no impact on our results, as we 
demonstrate in §4.2. The steady decline in the number of 
satellites at larger radii is genuine, and simply reflects the 
fact that the number density of satellites decreases with in- 
creasing halo-centric radius faster than (e.g.. Beers & 
Tonry 1986; Carlberg, Yee & EUingson 1997a; van der Marel 
et al. 2000; Lin, Mohr & Stanford 2004; van den Bosch et 
al. 2005; Chen 2007). Finally, there is also a clear indication 
that satellites at larger halo-centric radii are, on average, 
bluer. Again, this is consistent with numerous other studies 
(e.g., Biviano et al. 1996; CoUess & Dunn 1996; Carlberg et 
al. 1997b; Lares et al. 2004). 

The lower six panels of Fig. 2 show the average "'^ (p — r) 
colour of satellite galaxies as function of Mh (left column), 
Al, (middle column) and i?proj/i?i80 (right column). Lines 
of different colours correspond to different bins in one of 
the other two parameters, as indicated at the top of each 
panel. Errorbars (often barely visible) reflect the standard 
deviations obtained using the jackknife technique. To that 
extent we divide the group catalogue in TV = 20 subsamples 
of roughly equal size, and recalculate the mean colour 20 
times, each time leaving out one of the 20 subsamples. The 
jackknife estimate of the standard deviation then follows 
from 
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Figure 2. The contours in the upper panels reflect the l/Vmax-weighted distributions of satellite colours as functions of halo mass Mi^ 
(left), stellar mass M* (middle) and normalized, projected radius i?proj/^i80 (right). The connected solid dots indicate the average 
satellite colour as function of these quantities. The lower six panels show the average colour again as function of these three quantities, 
but this time split in bins of one of the other two quantities, as indicated. The values at the top of each panel indicate the central values 
of each bin. In particular, for -Rproj/^lSO we consider bins of [0.0,0.2] (green lines), [0.2,0.4] (blue lines), [0.4,0.6] (black lines), [0.6,0.8] 
(red lines) and [0.8,1.0] (magenta lines). For the stellar mass we adopt bins in log[A/»/(/i~^ 1^©)] of [9.0,9.5] (green lines), [9.5,10.0] 
(blue lines), [10.0,10.5] (black lines), [10.5,11.0] (red lines) and [11.0,11.5] (magenta lines). And flnally, for the halo mass we use bins 
in log[Mh/(h-l M0)] of [12.0,12.5] (green lines), [12.5,13.0] (blue lines), [13.0,13.5] (black lines), [13.5,14.0] (red lines) and [14.0,15.0] 
(magenta lines). Note that the last bin in is twice as broad as the other bins, which is done to assure a sufficient number of satellite 
galaxies to be able to calculate reliable statistics. Errorbars (sometimes barely visible) indicate the standard variance obtained from 20 
jackknife samples. Note how the average colours of satellite galaxies at flxed stellar mass are virtually independent of their environment 
(halo mass and halo-centric radius). 
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with Xi the average colour obtained from jackknife sample 
i. 

Panel (d) of Fig. 2 shows that at fixed stellar mass the 
average colour of satellite galaxies depends only very weakly 
on halo mass. On the other hand, at fixed halo mass, there is 
a very strong dependence on stellar mass, with more massive 



satellites being redder. The same trend is also visible from 
panel (h), which shows the average colour as a function of 
stellar mass for different bins in halo mass. Panels (f) and 
(g) show that the average satellite colour becomes slightly 
bluer with increasing halo-centric radius at fixed halo mass. 
At the same time, satellites at a fixed i?proj/i?i80 are also 
somewhat redder if they reside in more massive haloes. How- 
ever, as is evident from panel (e) , at fixed stellar mass there 
is no discernable dependence of the average satellite colour 

on iiproj/i?180. 
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Figure 3. Same as Fig. 2 but for tlie concentration parameter C = rgo/rso rather tlian the colour. Note that the overall trends are 
remarlcably similar. In particular, at fixed stellar mass the average concentration of satellite galaxies is virtually independent of their 
environment (halo mass and halo-centric radius). See text for a detailed discussion. 



Fig. 3 shows the same as Fig. 2 but now for the aver- 
age concentrations of the satellite galaxies, rather than their 
colours. Note that the overall behavior is remarkably similar 
to the case of the average colours: (i) at fixed halo mass there 
is a strong dependence on stellar mass (more massive satel- 
lites are more centrally concentrated), while at fixed stellar 
mass there is almost no halo mass dependence; (ii) at fixed 
halo mass there is a weak (barely significant) dependence on 
halo-centric radius (satellites at larger i?proj/^?i80 are less 
concentrated), and at fixed R-proj/Riso there is a similarly 
weak dependence on Mh; (iii) at fixed -Rproj/i?i80 satel- 
lite concentrations depend strongly on stellar mass, while 
at fixed stellar mass there is no discernable dependence on 

-Rproj/.Rl80. 

To summarize, the average colours and concentrations 
of satellite galaxies show almost no sign of intrinsic envi- 
ronment dependence. Rather, they seem to be completely 
determined by their stellar masses: to good accuracy Figs. 2 



and 3 can be used to read off the average colours and concen- 
trations of satellite galaxies of a given stellar mass mthout 
having to know anything about their environment. At least 
a significant part of the (already weak) dependence on halo 
mass is not causal, but merely owes to the fact that more 
massive haloes contain more massive satellites (cf. Fig. 1). 
Similarly, the mild dependence on Rproj/Riso seems to owe 
mainly to mass-segregation, rather than to a causal rela- 
tion between colour and halo-centric radius. This remark- 
able dearth of environment dependence for satellite galaxies 
is the main results of this paper. As we discuss in §6, it 
has profound implications for theories of galaxy formation. 
The remainder of this paper is therefore concerned with ad- 
dressing the robustness and statistical significance of these 
findings, and with reconciling them with previous, at first 
sight antithetical, results. 
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Figure 4. The standard deviations ap'^{g — r)] of the colour distribution of satellite galaxies as function of halo mass (left panels), 
stellar mass (middle panels) and projected radius (right panels). Lines of different colour correspond to different bins in one of the other 
three conditionals, as indicated. The values at the top of each panel indicate the central values of each bin, and the colour coding is the 
same as in Fig. 2. Errorbars reflect the standard deviations obtained from 20 jackknife samples. 



3.2 Higher order moments 

The above analysis, based on the average colours and con- 
centrations, suggests that the star formation histories and 
morphologies of satellite galaxies are virtually independent 
of their environment. Instead, they seem to be governed al- 
most entirely by their stellar masses. However, the averages 
only reflect the first moments of the full probability distribu- 
tions P{° \g - r)\M^,Mh,Rproj) and P{C \M,, Mh, Rproj). 
Since two distributions can be very different and yet have 
the same average, the above conclusion is only valid for the 
first moments of the colour and concentration distributions 
of satellite galaxies. 

As a logical next step, we therefore now focus on 
the second moments. Fig. 4 shows the standard deviations 
cl^'^ {g — r)] as functions of the three conditionals. The over- 
all behavior is very similar as for the averages: (i) at fixed 
halo mass there is a strong dependence on stellar mass, while 
at fixed stellar mass, the dependence on halo mass is much 
weaker, (ii) at fixed J?proj/^?i80 there is a strong depen- 
dence on M«, but at fixed stellar mass there is virtually no 
dependence on projected radius. Thus the second moments 
of PC'-'^ {g - r)\M, , Mh, Rproj) depend most strongly on M,, 
with only a weak dependence on M^, and virtually no de- 
pendence on Rproj/Riso- Fig. 5 shows the same but now 
for the satellite concentrations. Contrary to the colours, the 
standard deviations of P{C\Mt,, M^, Rpioj) seem to be al- 
most independent of stellar mass. In fact, within the errors. 



there is no really significant dependence on any of the three 
conditionals. 

We can therefore conclude that the statement that the 
star formation histories and morphologies of satellite galax- 
ies are (virtually) independent of their environment remains 
valid up to the second moments of their respective distribu- 
tion functions. 

Expanding this investigation, we could next present 
similar plots for the third moments, the fourth moments, 
etc. However, already for the third moment (or equivalently, 
the skewness) we find that the results become too noisy 
(i.e., the jackknife errors are too large) to make any quan- 
titative assessment. Therefore, we choose to show (a sub- 
set of) the full probability distributions instead. The up- 
per panels of Fig. 6 show the full colour distributions of 
satellite galaxies of a given stellar mass (indicated at the 
top of each panel) for three different bins in halo mass (in- 
dicated by the values in square brackets). The lower pan- 
els, on the other hand, show the colour distributions of 
satellites in haloes of a given mass for three different bins 
in stellar mass. The small vertical bars near the top of 
each panel indicate the averages of the corresponding dis- 
tributions. Only two distributions are shown in the upper 
right-hand and lower left-hand panels, since haloes with 
12 < log[Mh/{h~^ Mq)] < 12.5 do not contain satellite 
galaxies with log[M*/(/i~^ M©)] > 11. As is clear from a 
comparison of the upper and lower panels, the colour distri- 
butions of satellite galaxies are far more dependent on stellar 
mass than on halo mass. However, it is also clear that satel- 
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Figure 5. Same as Fig. 4 but for the concentrations of satellites rather than their colour. 



lite galaxies of the same stellar mass that reside in haloes 
of different masses have colour distributions that are not 
exactly the same. In fact, there are significant differences, 
though they typically relate to higher order moments of the 
distributions. For example, in the upper left panel, the distri- 
butions for \og[Mh/{h~^ M0)] > 13.0 are clearly negatively 
skewed, while that for 12 < log[Mh/(/i"^ M©)] < 12.5 has 
a positive skewness. We therefore conclude that the depen- 
dence of satellite colour on environment is extremely mild, 
and only apparent from the higher-order moments of the 
colour distributions. 

Finally, Fig. 7 shows the distributions of satellite con- 
centrations for the same bins in halo mass and stellar mass as 
in Fig. 6. The conclusions are the same as for the colours: at 
fixed stellar mass the distributions only depend very weakly 
on halo mass (even more weakly than in the case of the 
colour distributions) , while at fixed halo mass, the distribu- 
tions depend very strongly on stellar mass. 

4 TESTING THE ROBUSTNESS 

In order to check the robustness of the above results, we 
perform a number of tests, which we now describe in turn. 

4.1 Group masses 

As mentioned in §2 each group in our catalogue has two es- 
timates for the mass of its corresponding dark matter halo: 
one based on its ranking of the characteristic luminosity, 
defined as the total luminosity of all group members with 
"■^Mr — 5 log h < —19.5, and one based on its ranking of the 



characteristic stellar mass, defined as the total stellar mass 
of all group members with "'^Mr — 5 log h < —19.5 (see Y07 
for details). Thus far we have used the latter. Since the stel- 
lar masses are determined from the colours of the galaxies, 
one could potentially be worried that this has somehow in- 
fluenced the dependence of satellite colour on halo mass. To 
test this, we repeat our analysis using the halo masses based 
on the luminosity ranking. We also change the definition of a 
satellite galaxy: previously a satellite galaxy was defined as 
any group member that is not the most massive, while now 
a satellite galaxy is any group member that is not the most 
luminous. In a few cases this turns a satellite galaxy into 
a central galaxy, and vice versa. The left-hand column of 
Fig. 8 shows the results thus obtained, and which should be 
compared to the middle column of Fig. 2. Since we only con- 
sider satellite galaxies that reside in groups with an assigned 
halo mass in the range 12.0 < log[A/h/(ft-^M0)] < 15.0, the 
satellite samples in both cases are not exactly the same. In 
fact, when using the group masses based on the luminosity 
ranking we obtain a sample of 60835 satellites, compared to 
60245 in the case of the stellar mass ranking. This difference, 
however, is extremely small, so that the upper two panels 
in the first column of Fig. 8 are almost identical to those 
in the middle column of Fig. 2. The more interesting com- 
parison regards the lower panels in these columns, where 
the average colours as function of stellar mass are shown 
for different bins in halo mass. As can be seen, the results 
are almost indistinguishable, indicating that our results are 
extremely robust to whether we use the halo masses based 
on the luminosity ranking or the stellar mass ranking. We 
have verified that this is also the case for all other statistics 
presented in this paper. 
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4.2 Fiber Collisions 

As discussed in Y07, the galaxy sample on which group 
Sample II is based suffers from incompleteness due to fiber 
collisions: no two fibers on the same SDSS plate can be 
closer than 55 arcsec. Although this fiber collision constraint 
is partially alleviated by the fact that neighboring plates 
have overlap regions, ~ 7 percent of all galaxies eligible 
for spectroscopy do not have a measured redshift (Blan- 
ton et al. 2003d). Since fiber collisions are more frequent 
in regions of high (projected) density, they are more likely 
to occur in richer groups, thus causing a systematic bias. 
In Sample III we (partially) correct for this by assigning 
galaxies that lack an observed redshift due to fiber collisions 
(hereafter fiber-collision galaxies) the redshift of the galaxy 
with which it collided. As shown in Zehavi et al. (2002), 
roughly 60 percent of the fiber-collision galaxies have a red- 
shift within 500 km s^^ of their nearest neighbor, which jus- 
tifies the above procedure. However, there are also cases in 
which the fiber-collision galaxy has a true redshift that is 
very different from that of its nearest neighbor. In this case 
the above procedure will cause the group finder to falsely 
identify the fiber-collision galaxy as a group member (and 
with the wrong absolute magnitude and stellar mass) . In or- 
der to avoid these failures from impacting on our results, we 
have used Sample II as our fiducial sample in the present 
study (see §2). Nevertheless, since the above method seems 



to work well for ~ 60 percent of the fiber-collision galax- 
ies, it is worthwhile to check if our results are sensitive to 
whether we include fiber-collision galaxies or not. We there- 
fore repeat our analysis using Sample III. The selection cri- 
teria discussed in §2.1 now yield a sample of 92546 satellite 
galaxies (32301 more than for Sample II). Yet, the results, 
shown in the middle column of Fig. 8, are almost identical 
to those obtained from Sample II: once again, at fixed stel- 
lar mass the average satellite colour is almost independent 
of Mh and i?proj/^?i80- The main difference with respect to 
Sample II regards the tails of the colour-stellar mass distri- 
bution, as can be seen from the contour plot in the upper 
panel. In Sample III there is a tail of faint, very red galaxies, 
which is absent in Sample II. These galaxies are most likely 
fiber-collision galaxies for which the assigned redshift has 
a large error. Another difference, which is not shown here, 
is that in Sample III there is no decline in the number of 
satellites for i?proj/i?i80 0.15 as in the upper right panels 
of Figs. 2 and 3. However, since the fiber-collision galaxies 
do not have any preferred colour, this has no impact on the 
average colour of satellite galaxies at these small radii. 



4.3 Malmquist Bias 

As mentioned in §2.1, the sample of groups and galaxies used 
for our analysis is flux-limited, not volume-limited. In or- 
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Figure 7. Same as Fig. 6 but for the concentrations of satellites rather tlian their colour. 



der to correct for the consequential Malmquist bias we have 
weighted each galaxy by 1/Knax, the inverse of the comoving 
volume out to which the galaxy would have made it into our 
sample. However, in computing Vmax we did not account for 
the fact that the flux-limit of the SDSS varies mildly with the 
position on the sky (see for example Blanton et al. 2003b). 
In order to investigate how sensitive our results are to the 
1/Knax weighting scheme adopted, we now repeat our anal- 
ysis without any weighting. The results are shown in the 
right-hand column of Fig. 8. As expected the distribution of 
galaxies in the plane of colour vs. stellar mass is significantly 
different, showing a relatively larger contribution from mas- 
sive (and hence red) galaxies. As can be seen from the black 
solid line in the upper panel, at the low mass end satellites 
are about 0.1 magnitude bluer on average than in the case 
with weighting. A galaxy of a given stellar mass that is bluer 
will also typically be brighter, and can therefore be seen out 
to a higher redshift. In other words, by not applying the 
1/Knax weighting one artificially overestimates the fraction 
of blue galaxies at a given stellar mass, thus causing the av- 
erage colour to be too blue (but only mildly so). Yet, despite 
the fact that the average colours at a given stellar mass are 
somewhat different, even in this extreme case without any 
weighting to correct for Malmquist bias one finds that there 
is only a very mild halo mass dependence of (^'^{g — r)) at 
fixed M, and no significant dependence on R-proj/Riso at 
all, in good agreement with our fiducial results in Fig. 2. 

For brevity we only show the above tests for the average 
satellite colours as function of stellar mass. However, the av- 
erage colours as function of halo mass or halo-centric radius, 



as well as the average concentrations, all point to the same 
conclusion: the results presented in §3 are extremely robust. 
They do not depend on our choice of the group mass indica- 
tor, they are not infiuenced by sample incompleteness due 
to fiber-collisions, and they are not affected by an imperfect 
correction for Malmquist bias. 



5 COMPARISON WITH PREVIOUS STUDIES 

The analysis presented above shows that at fixed stellar mass 
the colours and concentrations of satellite galaxies are only 
weakly dependent on their environment. This differs sub- 
stantially from several previous studies, which have argued 
for a strong environment dependence of galaxy colours or 
star formation indicators (e.g., Hogg et al. 2004; Kauffmann 
et al. 2004; Blanton et al. 2005b; Weinmann et al. 2006a) 

Our study differs from most of these on the following 
grounds: 

(i) we have used halo mass and halo-centric radius as en- 
vironment indicators, as opposed to a projected number den- 
sity of galaxies. 

(ii) we have split our galaxy population in centrals and 
satellites and only focused on the latter. 

(iii) we have used broad-band colours, rather than more 
sophisticated star formation activity indicators, such as the 
equivalent width of various emission lines or the 4000Abreak 
strength, which are less sensitive to dust attenuation. 

The second of these issues is probably the most im- 
portant one. For example, as shown by Hogg et al. (2004) 
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and Kauffmann et al. (2004), the colour- magnitude relation 
of galaxies is strongly environment dependent, with lower 
density environments showing a far more pronounced blue 
sequence than high density environments. Our results, how- 
ever, suggest that the colour-magnitude relation of satellite 
galaxies has only a very weak environment dependence. We 
beUeve, though, that these results are not inconsistent with 
each other, but merely reflect that different environments 
have different satellite fractions: in high density environ- 
ments, such as clusters of galaxies, the (vast) majority of 
all galaxies are satellites. In low density environments, how- 
ever, most galaxies are central galaxies in low mass haloes. In 



paper I we have shown that central galaxies are, on average, 
bluer than satellite galaxies of the same stellar mass, and we 
believe that this is largely responsible for the environment 
dependence noticed in the studies mentioned above. 

The study whose results seems most antithetical to 
those presented here is that of Weinmann et al. (2006a; here- 
after W06), who used a SDSS galaxy group catalogue very 
similar to that used here in order to study the fractions of 
red and blue satellite (and central) galaxies as a function of 
halo mass and absolute r-band magnitude (see Hester 2006b 
for a similar analysis, though they did no split their sample 
in centrals and satellites). W06 find that at fixed absolute 
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magnitude the fractions of red and blue satellites are strong 
functions of halo mass, while at fixed halo mass there is only 
a very mild dependence on luminosity. At first sight these re- 
sults are completely opposite to the results presented here^. 
However, the following two reasons show that they are, in 
fact, consistent with each other. 

First of all, the fact that W06 found no significant lu- 
minosity dependence of the red fraction of satellites, /rod, 
owes entirely to the fact that they used a red/blue dividing 
line that is magnitude dependent. Clearly, for each magni- 
tude bin one can always identify a colour so that the satellite 
fraction red-ward of this colour is independent of luminos- 
ity. Apparently, this is roughly what happens when using the 
magnitude-dependent cut adopted by W06. liowever, since 
the dependence of colour on luminosity is now 'absorbed' 
by the magnitude dependence of the colour-cut, one can not 
use the absence of a luminosity dependence of /red to argue 
that colour docs not depend on luminosity. 

Secondly, the fact that W06 found a significant halo 
mass dependence for /red is consistent with our results. 
This is most easily seen from the upper left-hand cor- 
ner of Fig. 6. If we define the red satellite fraction, /rod, 
as the (l/Knax-wcighto(l) fraction of satellites red- ward of 
o i(gr-r) = 0.7, one obtains values of 0.33, 0.47, and 0.66 for 
the halo mass bins [12.0, 12.5] (red histogram), [13.0, 13.5] 
(black histogram) and [14.0, 14.5] (blue histogram), respec- 
tively. Thus, in agreement with W06, we find that /rod in- 
creases with increasing halo mass. Upon closer inspection, 
it is clear that this largely owes to the fact that the colour- 
distributions for these different halo mass bins have a dif- 
ferent skewness. However, as we have shown, the first and 
second moments are very similar. What really matters are 
the actual distribution functions, and there is no doubt from 
Fig. 6 that these depend more strongly on stellar mass than 
on halo mass. 

Wc therefore conclude that our results are not inconsis- 
tent with previous studies. We do caution, though, that one 
has to be careful not to overinterpret trends (or the absence 
thereof) in red or blue fractions. Since colour is an actual 
physical quantity, while /rod is not (and is sensitive to a 
somewhat arbitrary colour cut), we believe that our study 
presented here is more useful for gaining insight regarding 
the physical processes at work. 



6 CONCLUSIONS 

Using a largo SDSS galaxy group catalogue, we have exam- 
ined the ecology of satellite galaxies. In particular, we have 
studied how the colours and concentrations of satellite galax;- 
ies depend on their stellar mass and their environment. The 
latter is 'parameterized' in terms of the mass of the host 
halo and the halo-centric radius. Our main result is that 
both the colour and the concentration of a satellite galaxy 
are almost completely determined by their stellar mass, with 
only a very mild dependence on environment. Or, to put it 

t Although the results of W06 are based on r-band magnitudes, 
while those presented here are based on stellar masses, we have 
verified that our results hold if we were to use the r-band magni- 
tudes instead. 



in different words, in order to predict the (average) colour 
and concentration of a satellite galaxy, all one needs to know 
is its stellar mass. No information regarding its environment 
(halo mass or halo-centric radius) is required. 

There is a weak trend that satellites at smaller halo- 
centric radii and in more massive haloes are, on average, 
redder and more concentrated. However, we argue that nei- 
ther of these arc refiections of a causal environment depen- 
dence. Rather, they owe to mass segregation (more massive 
satellites reside on smaller halo-centric radii) and to the fact 
that more massive satellites on average reside in more mas- 
sive haloes. These two trends, combined with the fact that 
more massive satellites are redder and more concentrated, 
explain the weak environment dependencies evident in the 
data; they almost completely disappear when using stellar 
mass as a control variable (i.e., when only considering satel- 
lite galaxies in a narrow bin in stellar mass). 

It is interesting to combine these results with those of 
Paper I, in which we used the same SDSS group catalogue 
to investigate the colour and concentration differences be- 
tween central and satellite galaxies, in a statistical sense. 
That analysis has shown that satellite galaxies arc, on aver- 
age, redder and somewhat more concentrated than central 
galaxies of the same stellar mass. In addition, it weis shown 
that the average magnitude of this difference is independent 
of halo mass, in excellent agreement with the lack of en- 
vironment dependence presented here. Combining all these 
results, the following picture emerges: galaxies become red- 
der and somewhat more concentrated once they fall into a 
bigger halo (i.e., once they become a satellite galaxy). This is 
a clear manifestation of environment dependence. However, 
there is no indication that the magnitude of the transfor- 
mation (or its timoscale) depends on environment; a galaxy 
undergoes a transition when it becomes a satellite, but it 
does not matter whether it becomes a satellite of a small 
(Milky Way sized) halo, or of a massive cluster. In that re- 
spect there is no (significant) environment dependence, and 
this is the 'dearth' of environment dependence we are refer- 
ing to in the title of this paper. 

This dearth of environment dependence puts interesting 
constraints on the physical processes responsible for trans- 
forming satellite galaxies. As discussed at length in Paper I, 
the data are most consistent with a picture in which strangu- 
lation, i.e., the removal of the (hot) gas reservoir of satellite 
galaxies, is the main mechanism that operates on satellite 
galaxies, and that causes their transformation from the blue 
to the red sequence. The fact that satellites in clusters have 
the same properties as satellites of the same stellar mass 
in galaxy-sized haloes, strongly argues against mechanisms 
that are thought to operate only in very massive haloes, 
such as ram-pressure stripping or galaxy harassment. Note, 
though, that we are not claiming that these processes do not 
occur, after all, clear-cut example of ram-pressure stripping 
have been observed (e.g., Gavazzi et al. 1995; Kcrmey, van 
Gorkum & VoUmer 2004). We mainly argue that they are 
not the dominant processes causing satellites to undergo a 
blue-to-red sequence transition. We also emphasize that we 
have used the term 'ram-pressure stripping' to refer to a 
(rapid) stripping of the entire cold gas reservoir of a galaxy. 
Although we claim that this can not be the main cause 
for satellite transformations, our data is perfectly consistent 
with a picture in which ram-pressure stripping operates on 
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the cold gas, but only manages to strip the outer layers of 
the gas disk (i.e., the gas at large galacto-centric radii, where 
no or very little star formation occurs) . Such a stripping, af- 
ter all, has basically the same effect as strangulation (the 
removal of the hot gas) , as it effectively removes the fuel for 
future star formation (Hester 2006a,b). In fact, the observed 
HI deficiency of cluster galaxies (e.g., Warmels 1988; Cay- 
atte et al. 1990; Bravo-Alfaro et al. 2000; Levy et al. 2007) 
strongly suggests that a mechanism is stripping some, but 

not all, of the HI gas^. 

Our results also argue against 'pre-processing' in groups 
as an important process: as argued above, the only way in 
which a cluster is a 'special' environment, is that it has a 
larger fractional population of satellites than the 'field', and 
that the average satellite mass is higher than in the field. It 
has been suggested, though, that cluster galaxies are differ- 
ent from their field counterparts because they were part of 
an intermediate mass group prior to falling into the cluster 
(Zabludoff & Mulchaey 1998; Zabludoff 2002). Since merg- 
ing is more effective in groups than in either clusters or the 
field, this "pre-processing" could explain the morphology- 
density relation. However, our results show that there is no 
difference between satellite galaxies in clusters and satellite 
gala^xies of the same stellar mass in any other environment. 
This indicates that pre-processing in groups can not be the 
dominant cause for differentiating the cluster galaxy popu- 
lation from that of the field. A similar conclusion was re- 
cently reached by Berrier et al. (2008), who used numerical 
simulations to show that the vast majority of cluster galax- 
ies (~ 70%) fall into the cluster potential directly from the 
field, without having been a satellite in an intermediate mass 
group environment. 

We conclude by postulating, based on our results, that 
galaxy properties depend predominantly on stellar mass. 
The second most important 'parameter' for a galaxy is 
the distinction between centrals and satellites. When using 
these two parameters as control variables, there is virtu- 
ally no environment dependence left. This conclusion has 
recently received support from various directions. Kauff- 
mann et al. (2004) have shown that at fixed stellar mass, 
there is nearly no dependence of structural properties like 
Sersic index or concentration parameter on local galaxy den- 
sity. Pasquali et al. (2007) have shown that the trends of 
the isophotal shapes of elliptical galaxies ('disky' vs. 'boxy') 
with environment can be completely explained in terms of 
a pure stellar mass dependence. Using the same group cat- 
alogue as that used here, Pasquali et al. (2008) have shown 
that the star formation and AGN activity of galaxies has 
a much stronger dependence on stellar mass than on halo 
mass. Mouhcine, Baldry & Bamford (2007) find no depen- 
dence of the relationship between galaxy stellar mass and 
gas-phase oxygen abundance on local galaxy density. Blan- 
ton et al. (2005b) and KaufTmann et al. (2004) find that 
there is only a very weak dependence of galaxy size on local 



i Although this is often attributed to ram-pressurc stripping, vis- 
cous stripping (Nulsen 1997) or stripping due to tidal forces could 
also play an important role. In addition, the deficiency may also 
simply be a reflection of the fact that star-formation has used up 
a large fraction of the HI gas, without being replenished with new 
gas since its hot gas reservoir has been removed. 



density at fixed luminosity, which, as shown by Weinmann 
et al. (2008), is mainly a reflection of a difference between 
centrals and satellites. 
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